The phase diagram of the lattice Calogero-Sutherland model 
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We introduce a lattice version of the Calogero Sutherland 
model adapted to describe 1/d 2 pairwise interacting steps 
with discrete positions on a vicinal surface. The configura- 
tional free energy is obtained within a transfer matrix method. 
The full phase diagram for attractive and for repulsive inter- 
action is deduced. For attraction, critical temperatures of 
faceting transitions are found to depend on step density. 

PACS numbers: 68.35.Rh, 82.65. Dp, 68.35.Md, 64.60.Cn 



Among one dimensional quantum models the Calogero 
Sutherland (CS) model has attracted much attention 
In this model fermions on a line interact with a 
pairwise 1/d 2 potential. As an important application, the 
CS model was shown to be equivalent to a model of steps 
with 1/d 2 interaction, although in a continuous version, 
i.e. with no underlying lattice H. This approach pro- 
duced numerous useful results among which the analytic 
form of the free energy in the repulsive case (l],||,|[ . This 
exact result exerted a major influence on subsequent ap- 
proaches, nonetheless the CS model is a continuous one, 
with no crystallographic sites. As a consequence it can 
exhibit no roughening transition. In addition, the model 
is ill defined below a certain critical temperature To when 
the interaction is attractive ||. In this case, although 
there are mean field indications that faceting occurs at 
low temperature no solution for a better behaved 
model was suggested, nor were lattice effects taken into 
account when using the CS model in the framework of 
vicinal surfaces. 

In this article, in order to obtain a generic approach to 
study the instabilities of vicinal surfaces, we introduce a 
two dimensional model for steps which can be seen as a 
lattice CS model || . Since the lattice introduces a min- 
imal distance between steps (or fermions) this model is 
well behaved for all temperatures. For repulsive inter- 
action it unsurprisingly exhibits the Kosterlitz Thouless 
(KT) roughening transition. For attractive interaction 
we find first order faceting transitions which depend on 
step density. Using the transfer matrix method and finite 
size scaling we determine the full orientational phase di- 
agram in the repulsive and in the attractive case as well. 

In the field of surface physics the technological appli- 
cations of nanostructures make the understanding of the 
thermodynamical stability of surfaces a crucial issue. At 
a given temperature it is the behavior of the surface free 



energy as a function of orientation which determines the 
stability/instability of any surface orientation, as well as 
the equilibrium crystal shape (ECS) HQ. On the ex- 
perimental side the ECS was determined in several cases 
|^| [i0| . Orientational phase diagrams of Si, Pt and Au 
were drawn from x-ray experiments |]il|-|l3f . On the the- 
oretical side there have been several attempts to describe 
the energy of stepped surfaces in terms of pairwise inter- 
actions between steps. The form of step-step interactions 
is known in the case of elastic repulsion, which varies as 
A/d 2 jlj]. It is also known when electronic dipoles are lo- 
cated on steps, then the interaction is also 1/d 2 , although 
with either sign, according to dipole orientation. 

Although complex effects such as terrace reconstruc- 
tion are likely to play a role in some experimentally ob- 
served cases of faceting, by using the CS model, we focus 
on generic aspects of instabilities. 

Our model depends on two parameters, the kink cre- 
ation energy Ek and the interaction amplitude A. For 
N steps separated on average by Iq sites and with peri- 
odic boundary conditions the configurational energy for 
a periodic cell is written as : 
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where m and m' vary from 1 to N, h m ^ y is the position 
of step m at site y along the step direction. L = NIq is 
the size of the periodic cell. 

The first term of TL in Eq.(|l|) is the contribution of 
kinks. The next two terms reflect pairwise step step 
A/d 2 interactions. The potential V was introduced in 
order to include interactions between steps belonging to 
different periodic cells, in addition to the interaction be- 
tween steps belonging to a given cell |l) . The third term 
is the interaction of steps in the cell of interest with 
their periodic images, this term is absent in V. When 
taking the Boltzmann constant equal to unity the g pa- 
rameter of the CS model is g(T) = 2A/Tb 2 (T), where 
b 2 (T) = (cosh(E K /T) - l)- 1 is the T dependent step 
diffusivity. 

We determine here the orientational phase diagram of 
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the model i.e. critical temperatures in the (p, T) plane. 
The step density p = 1/Iq varies between and 1. In the 
following we take integer values for l Q . 

In order to localize transitions we perform a Finite Size 
Scaling (FSS) analysis of the free energy obtained with 
the Transfer Matrix (TM) method on a cylinder (see e.g. 
[[L5IQ). We thus treat N (with N from 2 to 6) infinitely 
long steps with transverse periodic boundary conditions, 
the free energy depend on the largest eigenvalue of the 
TM. Our goal in using FSS is twofold. On one hand, from 
the theory of conformal invariance, the variation of the 
free energy with N is a criterion for critical states |L7|. 
On the other hand, the behavior of the (N — > oo) extrap- 
olated free energy is a criterion for the thermodynamic 
stability of phases. For the FSS analysis of rough states 
it is of interest to introduce a coarse grained energy as a 
Fourier expansion on wave vectors q = (q m , q y ) Gq|: 
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where r\ m (resp. r\ y ) is the stiffness in the m (resp. y) 
direction. For this Gaussian model the FSS properties 
can be readily determined and the free energy per terrace 
site can be written : 



F{N) = F a 



(3) 



In the perspective of conformal field theory the central 
charge c is unity, since the system is rough in the sense of 
KT. The sound velocity v is thus directly related to the 
anisotropy in the m and y directions by v — pWVm/Vy 
(in agreement with p9[). 

Then, for a given p, we examine the convergence of 
the free energy obtained with the TM method by fitting 
numerical values with : 



F(N) = F a 



a 
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For g > 0, rough states are critical and a is expected 
to be equal to 2. Fig. (Q) shows the values of the appar- 
ent exponent a resulting from the fit for p = 0.25. In the 
region where a = 2 the system is rough. When using the 
coarse grained energy of Eq. (EJ) for a finite cylinder with 
approximate values of r\ m and r\ y |H],^(J to estimate v, 
we find for the apparent central charge c = 0.92 ± 0.02. 
With exact values of r\ m and r\ y the real central charge 
would thus be nothing but unity. In contrast, when the 
apparent value of a is greater than 2 the system is flat. 
The roughening temperature is located at the threshold. 
This can be confirmed by other types of arguments. Thus 
we have calculated with the TM method the standard de- 
viation of the terrace width distribution w 2 , which takes 
the special value 4j at the roughening temperature Tr 
p9| . Tr can also be determined by solving 



which is in principle valid for small p. It is apparent 
on Fig. (|j) that all criteria lead to very close values of 
Tr and we have made the phase diagram of Fig. (||) with 
the w 2 criterion. 
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FIG. 1. For p = 0.25 the free energy of the model of Eq. (0) 
for a finite cylinder is used to produce successive values of 
c(N), the plotted values of c are power law extrapolations. 
The values of a axe obtained by fitting data with Eq. (Q). 
Within the capillary approximation w 2 is equal to 4/ it 2 at 
Tr. 

Note that in the fermionic approach and subse- 
quent papers, the interaction of steps with their pe- 
riodic images in other cells was missing (third term 
in Eq. (|l])). This has led to two different determina- 
tions of the central charge, c = 1 was believed to be 
the exact value |2Ip^ ] while from Q one gets c — X 
with A = (I + v 7 ! + 2.g)/2. When the missing term 
= if H^IX^ 1S mc l u ded the fermion energy per 
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site becomes Ef = 

locity is v = 2 f^, this shows that the central charge is 
indeed unity. Neglecting this term would have made the 
FSS analysis inconsistent. 

As step-step interactions are pairwise one can invoke 
the step-hole symmetry of the Hamiltonian in Eq. (jl]) 
to show that the phase diagram exhibits a p — ► 1 — p 
symmetry. The criterion in Eq. (^) can thus be used for 
holes and the phase diagram is symmetric with respect 
to p = 1/2 . One sees that Eq. (||) gives an acceptable 
estimation for Tr even for p — 1/2, where it crosses the 
corresponding criterion for holes. Nonetheless it turns 
out that g(Tf() indeed varies when varying A/Ek, so 
that the scaling of Eq. (|1) is weakly broken by lattice 
effects around p — 1/2. 

In order for the free energy to reflect the step-hole sym- 
metry and exhibit more readily its convexity/concavity 
properties we use F s (p) = F(p)—pF(l) withF(l) = ^A. 
It is then convenient to introduce a symmetric Landau 
like expansion of F s . It turns out that excellent fits of 
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F s (p) are obtained with the following form : 

m = b + b 2 ( P - l -f + h( P - l -t + b e ( P ~ \f 
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with the constraints that <p(Q) = and that the slope 

2 

for p — is fl—^r-A, with (3 the free energy per unit length 
for an isolated step. The coefficients hi are temperature 
dependent. In the rough phase (low p) we find a good 
agreement between extrapolated values of the free energy 
F(p) and the classical analytical form [Q-Q: 



lip) = h+Pp + Sp 3 
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where /o, the surface energy of the reference plane, is 
here taken equal to zero and S(T) is given by 



S(T) = —Tb 2 
b 



(T)X 2 with A = i (l + ^TT2^j (8) 



Note that when Iq is not an integer the roughening 
temperature can be much lower than the value given by 
Eq. © ||. 

In the attractive case (g < 0) it is known that fermions 
collapse for g < —1/2. It was conjectured that a tricrit- 
ical point would appear and that faceting would occur 
if the potential divergency at the origin were cut off || . 
In our lattice model steps cannot collapse and the orien- 
tational phase diagram can be investigated. There are 
two parameters, the kink creation energy Ek and the at- 
traction strength A, or equivalently the temperature To 
defined by 3 (To) = —1/2. We have made calculations for 
Tq/Ek — 0.20, 0.60 and 1.00, in order to see the influence 
of this ratio. 
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FIG. 2. Orientational phase diagram. Repulsion (g > 0) : 
A/E K = 0.014 (o), A/E K = 0.140 (o). Vertical lines indicate 
the flat regions. Dotted lines represent the criterion for Tr 
as in Eq. (ISf) for steps and for holes. Attraction (<? < 0) : 
Tq/Ek = 0.20 (A), 0.60 (□) and 1.00 (*). Dashed lines are 
guides for the eye. 



We find that for T > To one has a — 2 for all tested 
values of p. The free energy is convex, the entropic repul- 
sion is dominant, so that all orientations are stable and 
rough. For T < T one finds that a increases above 2 
when T crosses aj) dependent temperature T c (p). Correl- 
atively the orientational form of the free energy changes 
and F s (p) becomes concave around p = and p = 1 (see 
Fig. H|). This indicates that the related orientations [0,p c ] 
and [1 — p c , 1] are unstable with respect to faceting. Of 
course, within the thermodynamical limit (N — > 00) TM 
calculations would give directly rise to phase separation 
and the free energy would converge like 1/N due to the 
presence of walls. Except at very low T, this is not the 
case and it is clear that concave regions are stabilized by 
finite size effects. The values of the corresponding free 
energy is presumably influenced by such effects whereas 
the tie line is not, very much alike convex regions. From 
the fitted p— 1/2 power expansion of F s (p) we determine 
the boundaries of the two-phase (faceted) p range. For 
Tq/Ek = 0.20 and T /E K = 0.60 the unstable range 
increases with decreasing T. Below a critical tempera- 
ture Ti , F s (p) becomes concave in the whole p range and 
faceting occurs between p = and p — 1. Fig.(|^) shows 
extrapolated values of F s (p) for various temperatures in 
the vicinity of Tq . 




FIG. 3. Symmetric free energy F s (p) (see text) as a func- 
tion of step density p, Tq/Ek = 0.60. TM results extrapo- 
lated with Eq. (^) (symbols), fit with expansion as in Eq. (^) 
( ), tie line for T/T = 0.848 ( ). T t /T = 0.838. 

It is noticeable that, for Tq/Ek = 1.00, Ti is very close 
to To. Only in the very narrow range of temperature 
[Ti,T ] does faceting occur close to p = and p = 1. 
Moreover, as b 2 turns out to be negative for Tj < T < T , 
an unstable p range appears around p = 1/2, as shown 
on Fig. (W). Results for the three cases are summarized 
on Fig. (g). 

As theoretical predictions and the analysis of experi- 
mental results rather use the expansion of F in powers 
of p as F{p) = fo + a\\p\ + a 2 p 2 + a 3 |p| 3 + ... it is of 
interest to show the coefficients deduced from our bi. 
It is noticeable that, the higher Tq/Ek, the narrower the 
p range where F s (p) is correctly given by the first four 
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terms of the p expansion. 
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FIG. 4. Same as Fig. @ for T /E K = 1.00. Tie lines for 
T/T = 0.996 (— - -) and ( ). T x /T = 0.989. 

For Tq/Ek = 0.20 this limited expansion is valid for 
p < 0.5. Fig. (|^) shows the values of 02 and 03 as func- 
tions of T. Above To the value of 03 is in agreement with 
S(T) as given by Eq. (^) whilst around To lattice effects 
bring about strong deviations. Thus 03 remains positive 
slightly below To, then it becomes negative, which indi- 
cates that the metastable phase computed with the TM 
method is uniform in this T range. The coefficient 02, 
which fluctuates around zero above T , becomes nega- 
tive below T), thus driving the instability and showing 
that the tricritical point is very close to Tq. The presence 
of a quadratic term for T < Tq is in contrast with mean 
field (Hartree Fock) approaches which give no such term 
j| . For Tq/Ek — 0.60 and 1.00 the full power expansion 
should be used for p > 0.1 in the whole T range. 
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FIG. 5. Coefficients of the power expansion of F a (p): ai, 
a;i and S(T) as in Eq. (^) as functions of T, dashed and dotted 
lines are guides for the eye, Tq/Ek = 0.20. 

We have thus shown that lattice effects alter the main 
features of the continuous CS model. In the repulsive 
case, while the previously known roughening transition 
is recovered, we show evidence for the symmetry related 
roughening of holes. In the attractive case the lattice 
CS model is well behaved and the critical temperature 



for first order faceting becomes p (step density) depen- 
dent. For vanishing p (and 1— p) the critical temperature 
tends to To, the limiting temperature in the continuous 
CS model. In the (p, T) phase diagram (0, Tq) and (1, To) 
appear as tricritical points. This phase diagram should 
be of interest for all lattice realizations of the CS model. 

We wish to thank V. Pasquier and R.A. Romer for 
useful discussions. 
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